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ABSTRACT 

We report on 7-ray observations of the Crab Pulsar and Nebula using 8 months of survey data with the 
Fermi Large Area Telescope (LAT). The high quality light curve obtained using the ephemeris provided by 
the Nancay and Jodrell Bank radio telescopes shows two main peaks stable in phase with energy. The first 
7-ray peak leads the radio main pulse by (281 ± 12 ± 21) [is, giving new constraints on the production 
site of non-thermal emission in pulsar magnetospheres. The first uncertainty is due to 7-ray statistics, and 
the second arises from the rotation parameters. The improved sensitivity and the unprecedented statistics 
afforded by the LAT enable precise measurement of the Crab Pulsar spectral parameters: cut-off energy at 
E c = (5.8 ± 0.5 ± 1.2) GeV, spectral index of T = (1.97 ± 0.02 ± 0.06) and integral photon flux 
above 100 MeV of (2.09 ± 0.03 ± 0.18) x 10~ 6 cm -2 s _1 . The first errors represent the statistical 
error on the fit parameters, while the second ones are the systematic uncertainties. Pulsed 7-ray photons 
are observed up to ~ 20 GeV which precludes emission near the stellar surface, below altitudes of around 
4 to 5 stellar radii in phase intervals encompassing the two main peaks. A detailed phase-resolved spectral 
analysis is also performed: the hardest emission from the Crab Pulsar comes from the bridge region between 
the two 7-ray peaks while the softest comes from the falling edge of the second peak. The spectrum of the 
nebula in the energy range 100 MeV - 300 GeV is well described by the sum of two power-laws of indices 
T sy nc = (3.99 ± 0.12 ± 0.08) and r /c = (1.64 ± 0.05 ± 0.07), corresponding to the falling edge of 
the synchrotron and the rising edge of the inverse Compton components, respectively. This latter, which links 
up naturally with the spectral data points of Cherenkov experiments, is well reproduced via inverse Compton 
scattering from standard Magnetohydrodynamics (MHD) nebula models, and does not require any additional 
radiation mechanism. 

Subject headings: gamma rays: observations - pulsars: individual (Crab, PSR J0534+2200) - ISM: supernova 
remnants - ISM: individual: Crab Nebula - Fermi 
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1. INTRODUCTION 

The Crab Nebula belon gs to the cla ss of filled-center su- 
pernova remnants (SNR) (Green 2006), i.e. without any de- 
tected shell component, and is well studied in almost all wave- 
length bands of the electromagnetic spectrum from the radio 
(10~ 5 eV) to very high energy 7-rays (nearly 10 14 eV). It is 
held to be the archetypical pulsar wind nebula, even though 
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its physical and spectral properties are unique. It is associ- 
ated with the supernova explosion rep orted by Chinese as- 
tronomers in 1054 AD. Several models (Kennel and Coroniti 
(fl984h . Ide Jager and Hardingl (fl99l . Ide Jager et alJ (119961) 
and references therein) describe the photon production pro- 
cesses taking place in this nebula. Synchrotron radiation from 
high energy electrons in the nebular magnetic field is respon- 
sible for the observed spectrum from radio to MeV, while 
inverse Compton (IC) scattering of the primary accelerated 
electrons off the synchrotron photons, far infrared and Cos- 
mic Microwave Background (CMB) produces high energy 7- 
rays. While these two mechanisms seem to provide a reason- 
able description of the overall non-thermal radiation of the 
Crab Nebula, one cannot exclude possible deviations from 
this simplified picture, and lAtovan and A haronian ( 1 19961) pro- 
posed that significant production of high energy 7-rays by 
bremsstrahlung radiation of relativistic electrons could take 
place in the Crab filaments. This paper, reporting the re- 
sults of a precise spectral analysis of the Crab Nebula between 
100 MeV and 300 GeV, adds new elements to this discussion. 

At the center of the nebula lies the Crab Pulsar, 
PSR J0534+2200, one of the most energetic known pulsars 
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(spin down power of E = 4.6 x 10 38 erg s _1 ), located at a 
distance of (2.0 ± 0.2) kpc. Estimation of its characteristic 
age using its rotation period (P = 33 ms) and derivative 
(P = 4.2 x 10 -13 s/s) yields an age of 1240 years, close to 
the observational value. 

The Energetic Gamma-Ray Experiment Telescope 
(EGRET), on orbit from 1991 to 2001 on board of the Comp- 
ton Gamma Ray Observatory (CGRO), reporte d the high 
energ y detection of the Crab Nebula and Pulsar (jNolan et al. 
fl99l Ide Jager et all [19961 iKuiper et al.1 |200T1) . A more 
detailed stu d y of t he 7-ray emission was then provided by 
iFierro etafl (1 1998b . presenting a complete phase-resolved 
spectral analysi s of the EGRET data, and more recently by 
iPellizzoni et alJ (120091) . describing the first AGILE timing 
results on 7-ray pulsars, including the Crab. Pulsations of 
the Crab Pulsar were reported above 25 GeV by the Major 
Atmospheric Gamma-ra y Imaging Ch erenkov Telescope 
(MAGIC) collaboration dAliu et al.ll2008l) . with a light curve 
consistent with the one measured by EGRET. 

Observations of the Crab Pulsar in high energy 7-rays can 
provide strong constraints on the locati on of the 7-ray emit- 
ting r egions: above the polar caps (Daughertv and Harding 
1996), in the intermediate models like the slot gap 
(Musl imov and Hardind 120041). o r far from the neutron star 
in the outer gaps (lRomanilll996h . In particular, the spectral 
analysis and the phase-resolved behaviour examined in this 
paper may be used to discriminate between these models. 

Successfully launched on June 11, 2008, the Large Area 
Telescope (LAT), aboard the Fermi Gamma-ray Space Tele- 
scope, formerly GLAST, offers the unique opportunity to 
study the high energy behaviour of the Crab Pulsar and Neb- 
ula in great detail. In this paper, we report the results of the 
analysis of the Crab region using 8 months of survey observa- 
tions with the Fermi-LAJL. In Sections|2]and[3] we describe the 
radio and 7-ray observations used, while Section [4] presents 
the results obtained from a detailed timing and spectral anal- 
ysis of the LAT data. Finally, in Sections [5] and [6] we dis- 
cuss and summarize the main implications of these results for 
models of both the pulsar and the nebula. 

2. RADIO TIMING OBSERVATIONS 

The Crab Pulsar, like many young pulsars, is affected by 
significant timing noise and glitches in rotation rate. Because 
of the long time interval considered in this article, the rota- 
tional phase behaviour with time has to be known with ex- 
treme precision. The Crab Pulsar is one of the more than two 
hundred pulsars of large spin-down luminosity E mon itored 
by the LAT pulsar timing campaign (S mith etai1l2008l) coor- 
dinated among Fermi, radio and X-ray telescopes. 

The timing solution for the Crab Pulsar has been built 
using ob servations made with the Nangay radio telescope 
(France) dTheureau et al.1 120051) and the J odrell Bank Obser- 
vatory 42 foot M KIA telescope (England) dHobbs et al.ll2004t 
iLyne et al.l[l993l) . A total of 698 times of arrival (TOAs) has 
been recorded between June 20, 2008 and April 8, 2009. The 
radio TOA dataset comprises 210 observations at 1.4 GHz 
and 488 observations at 600 MHz, in order to constrain the 
dispersion measure during the interval of the 7-ray observa- 
tions. This quantity is known to be highly variable, because of 
moving filaments in the Crab Nebula which make the column 
density of electrons along the lin e-of-sight change i n time. 

The TEMP02 timing package dHobbs et al.l2006l) was used 
to build the timing solution from the 698 TOAs. The mean 
time of arrival uncertainties are 2.66 and 15.64 [is for the 



Nangay and the Jodrell Bank Observatory observations, re- 
spectively. We fit the TOAs to the pulsar rotation frequency 
and its first two derivatives, as well as to the dispersion mea- 
sure (DM) and its first derivative (DM1) to take the variation 
of the electron column density into account. The fit further 
includes 10 harmonically-related sinusoids, using the "FIT- 
WAVES" option in the TEMP02 package, to flatten the tim- 
ing noise. We obtain DM = (56.7037 ± 0.0003) cm' 3 pc and 
DM1 = (3.05 ± 0.10) x 10~ 2 cm~ 3 pc yr' 1 on November 
11, 2008. The post-fit rms is 21.1 ^is, allowing for analyses of 
the 7-ray pulse profile with unprecedented precision. 

3. LAT DESCRIPTION AND OBSERVATIONS 

The LAT is an electron-positron pair conversion telescope, 
sensitive to 7-rays with energies from below 20 MeV to 
more than 300 GeV. It consists of a high-resolution con- 
verter tracker (direction measurement of the incident 7-rays), 
a CsI(Tl) crystal calorimeter (energy measurement) and an 
anticoincidence d etector to discrimin ate the background of 
charged particles (lAtwood et alJl200"9l) . In comparison to its 
predecessor EGRET, the LAT has a larger effective area (~ 
8000 cm 2 on-axis), a broader field of view (~ 2.4 sr) and 
a superior angular resolution (~ 0.6° 68% containment at 1 
GeV for events converting in the front section of the tracker). 

The following analysis was performed on 248 days of data 
taken in survey mode (August 02, 2008 - April 07, 2009). 
Events from the "Diffuse" class are selected, i.e. the highest 
quality photon data, having the most stringent background re- 
jection. In addition, we exclude the events with zenith angles 
greater than 105° due to the Earth's bright 7-ray albedo. 

4. RESULTS 

4.1. Light curves 

The selected 7-rays were phase-folded using the timing so- 
lution described in Section [2] Photons with an angle 9 < 
Max(6.68 — 1.76Logio(EM C v), 1-3)° of the radio pulsar po- 
sition, RA. = 83.63322°, Dec. = 22.01446° (J2000), are se- 
lected. This choice takes into account the instrument perfor- 
mance and maximizes the signal to noise ratio. At high ener- 
gies, the background is relatively faint compared to the Crab 
emission, so that a radius larger than the Point Spread Func- 
tion (PSF) can be kept. 

Using this energy-dependent region, the 7-ray light curve 
above 100 MeV is presented in Figures Q] and [2] (g). We 
have 22601 7-rays among which we estimate 14563 ± 240 
pulsed photons after background subtraction. Phase his- 
togra ms in radio (from the Nancay radio telescope) , opti- 
cal (lOosterbroek et alj|2008h . X-rays dRots et alJl2~0()4). hard 
X-ravs dMineo et al J 120061)7 7-ravs (CGRO COMPTEL and 
EGRET, Kuiper et al.l d2001 | )) an d very high energy (VHE) 7- 
rays (MAGIC. lAliu et alJ d2008l) ) are also plotted in Figure |2] 
We did not search for any correlation between giant pulses 
and 7-ray photons. 

The phase is taken at the maximum of the main radio peak 
observed at 1.4 GHz, as seen in Figures Q] and [2] (a). Consid- 
ering all events between 100 MeV and 300 GeV, two clear 
peaks PI and P2 can be seen at phase </>i = 0.9915 ± 0.0005 
and (j>2 = 0.3894 ± 0.0022, respectively. Hence, the peaks are 
separated by <5</> = 0.398 ± 0.003. PI and P2 are asymmet- 
ric. Their shapes can be well modeled by two half-Lorentzian 
functions (with different widths for the leading and trailing 
sides). The first peak presents rising and falling edges of half- 
widths 0.045 ± 0.002 and 0.023 ± 0.001 respectively. P2 
shows a slow rise and a steeper fall. The rising and falling 
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Phase 

FIG. 1 . — Light curve obtained with photons above 100 MeV within an energy-dependent circular region, as described in Section l4~71 The light curve profile is 
binned to 0.01 of pulsar phase. Insets show the pulse shapes near the peaks, binned to 0.002 in phase. The radio light curve (red line) is overlaid (arbitrary units). 
The main peak of the radio pulse seen at 1 .4 GHz is at phase 0. Two cycles are shown. 
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edges of P2 have Lorentzian half- widths of0.115±0.015 and 
0.045 ± 0.008 respectively. Hence, the 7-ray first peak leads 
the radio main pulse by phase 0.0085 ± 0.0005, as shown 
in Figure [T] where the radio profile (red line) is overlaid for 
comparison. 

The second 7-ray peak leads the second 1.4 GHz radio 
pulse (interpulse) by 0.0143 ± 0.0022 in phase. The peak 
separation is slightly wider at 1 .4 GHz than in 7-rays. 

An error in these 7-radio delays can also arise from the 
measureme nt of the dispersion measur e and its derivative. 
Following Man chester and Taylorl(ll977l) . the error on the dis- 
persion delay in the propagation of a signal at a frequency / 
through the interstellar medium is: 



A(At) 



ADM 



(1) 



where ADM takes into account the error on the measure- 
ment of DM and its derivative, and K = 2.410 x 10~ 4 
MHz -2 cm -3 pc s _1 is the dispersion constant. This yields 
a formal uncertainty of 1 .4 ^s, which is significantly smaller 
than the 21.1 /is accuracy of the overall timing solution, and 
therefore leads to an error of 0.0006 in phase on the 7-radio 
delay. 

The presence of a radio feature referred to as Low Fre- 



quency Component (LFC) by iMoffett and Hankinsl (11996b 
can be noticed, at phase 0. 896 ± 0.00 1 on the radio light curve 
obtained at 1.4 GHz as seen in Figures [T]and|2](a). This peak 
is assumed to be near the closest approach of the magnetic 
axis. The first 7-ray peak lags the LFC by 0.095 ± 0.002 in 
phase. 

Figure |3] shows the light curves in 5 energy bands, covering 
the 100 MeV - 300 GeV interval while Table Q] reports the 
evolution of the positions of the peak maxima (<fii and <p2 for 
PI and P2 respectively) and their half- widths (HW), for the 
energy bins between 100 MeV and 10 GeV. The photon num- 
ber counts above 10 GeV were not sufficient to fit the peak 
profiles. The phases of the first (PI) and second (P2) peaks 
do not show any significant shift with energy. Both become 
narrower when the energy increases, showing in particular a 
steepening in the P2 falling edge. 

Table Q] also presents the energy dependence of the rela- 
tive weight of the two peaks. The diffuse and nebular back- 
ground photon density has been first estimated in the 0.52 - 
0.87 phase interval, then renormalized and subtracted so as to 
determine the number of pulsed photons in both peaks. PI and 
P2 are here defined in the 0.87 - 1.07 and 0.27 - 0.47 phase 
interva ls respectively. As for the Vela pulsar ( Ab do et al.l 
2009a), the ratio P1/P2 decreases with increasing energy, es- 
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FIG. 3. — Fermi light curves for the Crab Pulsar in different energy bands within an energy-dependent circular region, as described in Section |4~T1 The light 
curve profile is binned to 0.01 of pulsar phase, except above 10 GeV, which is binned to 0.02 in phase. Two cycles are shown. 
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pecially above a few GeV. 

We define the off-pulse window as the 0.52 - 0.87 phase 
range, due to the bright emission of the pulsar in the rest 
of the phase. In the light curve above 10 GeV, we can no- 
tice an enhancement indicating a potential third peak at phase 
~ 0.74, coincident with the radio peak observed between 4.7 
and 8.4 GHz and referred to as High Fr equency Component 
2 (HFC2) bv lMoffett and Hankinsl (119961) . The excess above 
the background level (estimated at 2.10 counts per bin, with 
a bin width of 0.02 in phase) is 13.8 photons in the off-pulse 
interval. The statistical significance of this third peak, 2.3 er, 
is therefore too low to claim a definite detection and a third 
peak will not be considered separately in the analysis of the 
Crab Nebula. 

Figure|4]shows the counts maps of pulsed and nebular emis- 
sion in a 15° x 15° region centered on the pulsar radio posi- 
tion, for different energy bands. The nebular emission seen in 
the off -pulse window has been renormalized to the total phase 
(bottom row) and subtracted from the whole phase emission, 
to obtain the maps presenting the pulsed emission only (top 
row). The positions of the pulsar and nebula are coincident 
to within our angular resolution and the nebula appears as a 
point-like source. While the pulsed emission dominates in 
the on-pulse window, the nebula stands out in the off-pulse 
interval from the emission of the diffuse background at high 
energy only. 

4.2. Spectral analysis of the Crab Nebula 

The spectral analysis of the 7-ray emission of the Crab 
Nebula was perform ed using a maximum-likelihood method 
(Matt oxet al.lll996b implemented in the Fermi Science Sup- 
port Center science tools as the "gtlike" code. This fits a 
source model to the data along with models for the instru- 
mental, extragalactic and Galactic backgrounds. We used an 
updated instrument response function, Pass6_v3, that corrects 
a pileup effect identified in orbit. We selected photons in 
the 0.52 - 0.87 pulse phase window in a 20° region around 
the pulsar radio position. Owing to uncertainties in the in- 
strument performance still under investigation at low ener- 
gies, only events in the 100 MeV - 300 GeV energy band 
are analysed . The Galactic diffu se emission is modeled using 
GALPROP (IStrong et al.ll2004allbh updated to include recent 
HI and CO surveys, more accurate decomposition into Galac- 
tocentric rings, and many other improvements. The GAL- 
PROP run designation for our model is 54_59Xvarh7S. The 
instrumental background and the extragalactic radiation are 
described by a single isotropic component with a power-law 
shape. Sources nearby the Crab with a statistical significance 
larger than 5 a are ext r acted u sing the analysis procedure de- 
scribed in Abd o et ail (12009b!) but with 6 months of survey 
data, and taken into account in the study. 

The systematic errors on the spectral parameters are dom- 
inated by the uncertainties in the LAT instrument response 
functions (IRFs). We bracket the energy-dependent effective 
area with envelopes above and below the nominal curves by 
linearly connecting differences of (10%, 5%, 20%) at log(E) 
of (2, 2.75, 4) respectively. This yields the systematic errors 
cited below. 

In parallel to the standard analysis, we have also evaluated 
the sp ectrum using an unfolding method based on Bayes' the- 
orem (iD'AgostinilliggSHMazziottal^OOgl) . that allows the re- 
construction of the true energy spectrum from the observed 
one taking into account the dispersions introduced by the in- 
strument response function and without assuming any model 



for the spectral shape. The results from this analysis are con- 
sistent with those from the like lihood analysis . 

Using EGRET observations, de Jager et al. ( 1996) reported 
that the inverse Compton component dominates above ~ 200 
MeV whereas the synchrotron component is more significant 
at lower energies. Hence, the selected 7-ray photons should 
allow the study of both the fall of the synchrotron and the rise 
of the IC radiation. 

The spectrum of the Crab Nebula between 100 MeV and 
300 GeV is well described by the sum of two power-law spec- 
tra. As seen on the spectral energy distribution in Figure |5j 
one of the components decreases while the second increases 
with energy. We identify them as the falling edge of the syn- 
chrotron component and the rising edge of the IC component, 
respectively. The nebular spectrum can be modeled with the 
following function: 

^ = N sync {E GeV )- T ^ +N IC (E GeV )~ r ' c 

cnT 2 s^MeV -1 (2) 

where N sync = (9.1 ±2.1 ±0.7) x 10- 13 cm- 2 s" 1 MeV"\ 
N IC = (6.4 ± 0.7 ± 0.1) x 10~ 12 cm" 2 s" 1 MeV" 1 
are the prefactors determined on 35% of the total phase, 
T sy nc = (3.99 ± 0.12 ± 0.08) and T IC = (1.64 ± 0.05 ± 0.07) 
the spectral indices of the synchrotron and IC components. 
While the power-law index Tjc for the inverse Compton com- 
ponent provides a measure of the index of the mean elec- 
tron/positron energy spectrum in the nebula, the synchrotron 
index T sync possesses much less physical information, being 
just an indication of the steepness of the quasi-exponential 
turnover of the synchrotron component that peaks below the 
LAT energy window. Adopting a power-law fit to the syn- 
chrotron contribution apparent in the 100 - 400 MeV range 
is therefore a useful mathematical convenience. The corre- 
sponding flux above 100 MeV and renormalized to the total 
phase is (9.8 ± 0.7 ± 1.0) x 10~ 7 cm~ 2 s -1 . The first error 
is statistical, whereas the second is systematic. 

Figure shows the spectral energy distribution in E 2 ^ of 
the Crab Nebula renormalized to the total phase. The Fermi- 
LAT spectral points were obtained by dividing the 100 MeV 
- 300 GeV range into logarithmically-spaced energy bins and 
performing a maximum likelihood spectral analysis in each 
interval, assuming a power-law shape for the source. Above 
5.5 GeV, the width of the energy intervals is multiplied by 3 
to reduce the statistical uncertainties. These points, provid- 
ing a model-independent maximum likelihood spectrum, are 
overlaid with the fitted model described above over the total 
energy range (black curve). The fit of the synchrotron (purple 
dashed line) and IC (blue dashed-dotted line) are also repre- 
sented. The spectral points and the model agree well. Statis- 
tical errors (black error bars) and the overall error (red error 
bars) are plotted for the Fermi points. The EGRET spectral 
points are represented on th e same plot. As in the case of the 
spectrum of the Vela pulsar dAbdo et al.l2009al) . derived using 
an earlier set of response functions, Pass6_vl, markedly dif- 
ferent from Pass6_v3 at low energies, the LAT spectral points 
at high energy indicate a lower flux in comparison to EGRET. 
However, it can be noticed that the Fermi flux is higher than 
the EGRET flux, in the low energy band dominated by syn- 
c hrotron radiatio n. 

Ide Jager et alj (11996b found evidence in the EGRET data 
that the Crab synchrotron cut-off energy varied on time scales 
of the order of a year. We do not see significant variation 
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FIG. 4. — Counts maps (arbitrary units) presenting the pulsed (top row) and nebular (bottom row) emission, in three energy bands. Each panel spans 15° X 
15° in equatorial coordinates and is centered on the pulsar radio position. Left: 100 MeV < E < 300 MeV; Middle: 300 MeV < E < 1 GeV; Right: E > 1 GeV. 



TABLE 1 

Detailed parameters of the Crab Pulsar light curve. 



Energy interval 


<h 






4>2 






P1/P2 ratio 


(GeV) 


(xl0~ 2 ) 


(xlO" 2 ) 


(xl0~ 2 ) 


(xl0~ 2 ) 


(xl0~ 2 ) 


(xlO" 2 ) 




0.1-300 


99.2±0.1 c 


4.5 ± 0.2 C 


2.3 ± 0.1 c 


38.9 ± 0.2 C 


11.5 ± 1.5 C 


4.5 ± 0.7 C 


1.60 ±0.06 


0.1-0.3 


99.2 ±0.1 


6.0 ± 0.4 


2.3 ± 0.2 


38.3 ±0.8 


8.4 ±2.2 


8.1 ± 3.7 


1.73 ± 0.12 


0.3-1.0 


99.1 ±0.1 


4.3 ±0.2 


2.7 ± 0.2 


39.3 ± 0.3 


13.6 ± 1.7 


3.6 ±0.9 


1.60 ±0.08 


1.0-3.0 


99.2 ±0.1 


3.5 ±0.3 


2.3 ± 0.3 


38.2 ± 0.5 


8.5 ± 2.8 


6.0 ± 1.9 


1.49 ±0.12 


3.0-10.0 


99.5 ± 0.2 


2.6 ±0.5 


1.6 ±0.6 


38.9 ± 0.6 


5.3 ± 1.9 


2.0 ± 1.4 


0.95 ± 0.20 



a, b: These half-width (HW) parameters were obtained considering two half-Lorentzian distributions, for the rising and 
falling edges of the peaks respectively. 

c: These parameters were derived from a light curve binned to 0.002 of pulsar phase. 



in either the synchrotron or inverse Compton components in 
our more limited data span on time scales of one, two, or 
four months. As shown in Figure [3J a difference in flux 
is observed between EGRET and Fcr/wi-LAT in the en- 
ergy band dominated by synchrotron radiation as well as 
at higher energies (above 1 GeV). Even if variability in 
the synchrotron tail could be expected between EGRET 
and LAT, the lifetimes of the electrons producing gamma- 
rays via inverse Compton scattering are comparable to the 
remnant age, implying that the IC component should be 
steady in time. For these reasons, the flux change seen in 
the synchrotron component between EGRET and Fermi- 
LAT cannot be considered as significant. 

The photon counts at high energy are too few for a signifi- 
cant cut-off or break to be seen in the flux distribution of the 



IC component. No cut-off or break energy can be determined 
at low energy for the synchrotron component using the LAT 
data only. 

4.3. Spectral analysis of the pulsed emission 

Photons from both on- and off -pulse intervals are now con- 
sidered to analyze the pulsed emission. The spectral param- 
eters of the Crab Nebula mentioned in the previous section 
have been renormalized to match the total phase interval and 
fixed to perform the spectral analysis of the Crab Pulsar. 

After testing different functional forms to describe the spec- 
trum of the pulsar, we found the best fit to be given by an 
exponential cut-off power-law shape: 

^ = N (E GeV )- r e~ E / E < cm~ 2 s^MeV^ 1 (3) 
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where N Q = (2.36 ± 0.06 ± 0.15) x 10" 10 cm- 2 s" 1 MeV" 1 
is the prefactor, V = (1.97 ± 0.02 ± 0.06) the spectral index 
and E c = (5.8 ± 0.5 ± 1.2) GeV the cut-off energy of 
the distribution. The integral flux above 100 MeV is equal to 
(2.09 ± 0.03 ±0.18) x 10~ 6 cm~ 2 s" 1 . These results are 
consistent with the pulsed spectrum derived from the unfold- 
ing analysis. 

Figure [6] shows the spectral energy distribution of the Crab 
Pulsar over the whole pulse period compared to EGRET spec- 
tral points. Results of both experiments agree well in the 
100 MeV - 8 GeV energy range, i.e. even at low energies, 
where such consistency is not observed in the spectrum of the 
Crab Nebula. The larger energy band covered by the LAT and 
its better sensitivity allows us to determine the cut-off energy 
of the spectrum, which was not possible with EGRET. 

We also attempted to fit the data using a power-law with 

a generalized cut-off of the form er^ E l E ^ and found b = 
(0.89 ± 0.12 ± 0.28) with a likelihood value which is not 
significantly better than that obtained in the case of a simple 
exponential b = 1 cut-off. We compute the probability of in- 
correct rejection of other spectral shapes using the likelihood 
ratio test. For instance, if only statistical errors are included, 
the power-law and hyper-exponential b = 2 hypothesis shapes 
are rejected at a level of 10.7 a and 4.9 a respectively. 

4.4. Phase-resolved spectral analysis of the Crab Pulsar 

The large number of photons detected from the Crab allows 
a detailed phase-resolved spectroscopic study of its emission. 
Therefore, the pulse profile is divided in several intervals. The 
phase bins are chosen so as to contain ~ 1000 pulsed pho- 
tons in the energy-dependent region defined in Section 14.11 
A maximum-likelihood spectral analysis is performed in each 
pulse phase interval, assuming a power-law and an exponen- 
tial cut-off power-law shape to describe the pulsed emission. 

The definition of the phase intervals is given in Table |2] 
along with the spectral results. The last column lists the 
significance of the improvement obtained when using an ex- 
ponential cut-off power-law instead of a pure power-law, in 
terms of % 2 , if only statistical errors are included. 

The corresponding spectral energy distributions are pre- 
sented in Figure|7] where the horizontal error bars delimit the 
energy intervals. 90% C.L. upper limits were computed when 
the statistical significance of the energy interval was lower 
than 3 a. Figure [8] summarizes the phase-dependence of the 
variation of the spectral parameters, spectral index, cut-off en- 
ergy and integral flux above 100 MeV. The vertical error bars 
take into account both the statistical and systematic uncertain- 
ties in the spectral parameters, while the horizontal error bars 
delimit the phase intervals. 

One observes a slight steepening of the spectrum with phase 
in the interval corresponding to PI (0.955 - 0.098), with aver- 
aged values of spectral index and cut-off energy close to ~ 1 .9 
and ~ 3 GeV respectively. 

The 0.098 - 0.286 pulse phase interval presents the hardest 
spectrum with a spectral index of 1 .49 ± 0.09 ± 0.05. This re- 
su lt is consistent with the spectrum of the "bridge", as defined 
in lFierroetal1(fl998l) . 

The spectral indices of the two peaks are the same, within 
the error bars, but the second peak (0.286 - 0.410) is charac- 
terized by a cut-off energy apparently larger than that of PI. 
This difference is consistent with the decrease of the P1/P2 
ratio with the energy, especially above a few GeV, observed 
in Figure [3] 



Finally, the 0.410 - 0.520 phase bin has the softest spec- 
trum of the total pulse phase interval with a spectral index of 
2.28 ± 0.08 ± 0.10. This explains the trend seen in the phase 
histograms: in Figure|3] the right edge of the second peak falls 
with increasing energy. 

5. DISCUSSION 

5.1. Synchrotron and inverse Compton emission from the 
Crab Nebula 

The Crab Nebula is detected across the whole electromag- 
netic spectrum from radio to very high energy 7-rays. The 
total spectral energy distribution of this source is shown in 
Figure [9] from soft to very-high energy 7-rays. The spec- 
tral points obtained with the LAT data analysis are also repre- 
sented (red points). 

With a spectral index of T IC = (1.64 ± 0.05 ± 0.07), the 
LAT results on the rising edge of the inverse Compton com- 
ponent are consistent with EGRET (1.85+q f 5 , de Jag er et aTl 
(1996)). As can be observed in Figure [9] the highest part of 
the LAT spectrum links up satisfactorily to the lower energy 
bound of the Cherenkov data points. Using the LAT spectral 
parameters scaled to the full pulse phase, we obtain a flux at 
77 GeV of (1.18 ± 0.22 ± 0.37) x 10~ 14 cm" 2 s" 1 MeV" 1 
which agrees with the MAGIC differential flux at this energy 
of (1.14 ± 0.27 ± 0.34) x 10~ 14 cm" 2 s" 1 MeV" 1 . 
The Cherenkov and Fermi-LPC£ data now cover the entire in- 
verse Compton peak, as can be seen in Figure [9] and a break 
is expected at ~ 100 GeV. Although no significant cut-off is 
observed in the LAT data with the current statistics, the de- 
termination of its position with an increased Fermi-hPX data 
sample would help the calib r ation of Cherenkov telescopes, 
as discussed in lBastieri et alj (120051) . 

The inverse Compton scattering of relativistic electrons on 
the synchrotron, far infrared, and cosmic microwave back- 
ground radiation fields is considered to be the most prob- 
able mechanism for production of 7-rays above 1 GeV. 
However, using a sophisticated approach carried ou t in th e 
framework of the MHD flow of iKennel and CoronitH dl984l) . 
lAtoyan and Aharonianl ([1996) have commented on the appar- 
ent deficit of GeV photons in their calculations. Taking into 
account both EGRET and Cherenkov results and assuming 
a mean magnetic field which reproduces the very high en- 
ergy spectrum, they proposed that the high 7-ray flux ob- 
served by EGRET in comparison to their model is due to 
the enhancement of the bremsstrahlung emission from elec- 
trons captured in dense filaments. Figure [9] presents the 
broad-band energy spectrum of the Crab Nebula together with 
the inverse Compton model predictions from Figure 14 of 
lAtovan and Aharonianl (TT996) for three different values of the 
mean magnetic field for the nebula. In view of the results ob- 
tained with the LAT, modeling the data does not require any 
additional emission component. The Fermi-LAZI, in combi- 
nation with the Cherenkov observations above 100 GeV, are 
in good agreement with the 7-ray flux predicted from simple 
IC scattering when the magnetic field lies between 100 /iG 
and 200 /iG, i.e. below the canonical equipartition field of the 
Crab Nebula of 300 pG. This result is consistent with the es- 
timate of the magnetic field s trength B ~ 140 /iG obtained by 
iHorns and Aharonian (2004]). 

Concerning the low energy part of the nebular spectrum, 
the LAT spectral points, combined with COMPTEL's (tak- 
ing into account statistical errors only for the latter), can be 
fitted with a power-law with an exponential cut-off, follow- 
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TABLE 2 

Phase interval definitions and corresponding spectral parameters . 



1 


Wmax 


rlux 


Spectral index 


Cut-off energy (GeV) 


Ji.h 


0.870 


0.955 


21.0 ± 1.1 ± 3.3 


i ao i n i n i /~\ t\ 

2.03 ± 0. 12 ± 0.20 


1 ^7 I f\ C I A 1 

1.7 ± 0.5 ± 0.3 


5.1(7 


0.955 


0.971 


106.1 ± 6. 1 ± Is. I 


2.05 ± 0.07 ± 0.22 


^ C 1 A f 1 A f 

2.5 ± 0.6 ± 0.6 


6.0(7 


0.971 


0.981 


1 *"7"7 T 1 (1 O 1 1 1 A 

177. j ± 9.2 ± 14. U 


1 A"7 1 A m 1 A A"7 

l.y7 ± U.U/ zb U.U7 


TO 1 A £. 1 AO 

2.5 ± 0.6 ± 0.2 


7.2(7 


0.981 


0.987 


TTIO 1 AAA 1 1 c n 

232.8 ± 14.4 ± 15.0 


1 f\ A I AAO I A AC 

1.94 ± 0.08 ± 0.05 


2.7 ± 0.7 ± 0.2 


6.4cr 


0.987 


0.993 


264.0 ± 11.5 ± 11.7 


1.93 ± 0.06 ± 0.04 


4.3 ± 1.0 ± 0.4 


6.9ct 


0.993 


1.000 


205.0 ±7.5 ±27.1 


1.90 ±0.05 ±0.11 


5.5 ± 1.3 ±2.0 


6.6(7 


0.000 


0.016 


94.8 ±3.9 ±5.0 


1.84 ±0.08 ±0.03 


3.1 ±0.9 ±0.2 


6.9(7 


0.016 


0.098 


15.3 ±0.9 ±2.7 


1.74 ±0.07 ±0.20 


6.3 ± 1.8 ±2.5 


5.9cr 


0.098 


0.286 


5.8 ± 0.4 ± 0.4 


1.49 ±0.09 ±0.05 


5.5 ± 1.3 ± 0.7 


7.9cr 


0.286 


0.338 


25.5 ± 1.3 ±3.1 


1.72 ±0.08 ±0.14 


3.6 ± 0.8 ± 0.9 


6.8a 


0.338 


0.366 


52.3 ±2.2 ±2.5 


1.94 ±0.06 ±0.04 


6.2 ± 1.8 ±0.3 


5.0(7 


0.366 


0.386 


70.4 ± 2.8 ± 8.0 


1.92 ±0.05 ±0.04 


6.8 ± 1.8 ± 0.7 


5.8a 


0.386 


0.410 


65.8 ± 2.7 ± 13.3 


2.04 ±0.06 ±0.16 


10.0 ± 4.8 ± 11.6 


3.3cr 


0.410 


0.520 


13.7 ±0.8 ±1.8 


2.28 ±0.08 ±0.10 


7.3 ± 4.8 ± 2.4 


2.3(7 



a : in units of 10~ Y photons cm -2 s — 1 and divided by the width of the phase interval. 

b : this value characterizes, for each phase interval, the improvement obtained by using an exponential 

cut-off power-lawinstead of a pure power-law shape, to describe the pulsed spectrum. 



TABLE 3 

The radio delay with respect to other frequencies. 



Spectral band 


Radio delay 


Reference 


Optical 


255 ± 21 


(1) 


X-rays 


344 ± 40 


(2) 


Hard X-rays 


280 ± 40 


(3) 


7-rays (EGRET) 


241 ± 29 


(3) 


7-rays (LAT) 


281 ± 12 ±21 


(4) 



References: (1): Oosterbroek et al. (2008); (2): Rots et al. (2004); 



(3): Kuiper et al. (2003); (4): this paper. 

ing de Jag er et alJ {1996). The cut-off energy is estimated at 
E C}S ync = (97 ± 1 2) MeV. The h i gher value of this energy 
compared to that of Ide Jager et al.1 (1996) is due to the larger 
flux obtained with Fermi than by EGRET for the synchrotron 
component. The fit is represented with a blue dashed curve in 
Figure [9] 

5.2. High energy emission from the Crab Pulsar 

The high-quality statistics obtained with the Fermi-LAI 
both on the light curve and the spectrum of the Crab Pul- 
sar, allow a more detailed comparison with theoretical mod- 
els than previously possible. Currently, there are two classes 
of models that differ in the location of the emission re- 
gion. The first comprises polar cap (PC) models which 
plac e the emission near the magn etic poles of the neutron 
star dDaughertv and Hardin dfl996l) . The second class con- 
sists of the outer gap (OG) models dRom ani 1996), in which 
the emission extends between the null charge surface and 
the light cylinder, and the two-pole caustic (TPC) models 
(iDyks and R udak 2003|) w hich might be realized in slot gap 
(SG) acceleration models (Muslim ov and Ha rding 2004]), in 
which the emission takes place between the neutron star sur- 
face and the light cylinder along the last open field lines. 

Observations of the time delay between emission at differ- 
ent wavelengths have been reported previously: Table [3] sum- 
marizes the delay of the radio main pulse with respect to the 
first peak seen from optical to high energ y 7-rays. The LAT 
has a timing accuracy better than 1 /is dAbdo et al.l F2009c) 
and thus enables an accurate estimation of the abs olute po- 
sitions of the 7-ray peaks: it was shown in Section |4~T1 that 



the first 7-ray peak leads the radio main pulse by 0.0085 
± 0.0005 ± 0.0006 in phase, or (281 ± 12 ± 21) /is in 
time. Taking into account the presence of the LFC at phase 
0.896 ± 0.002, the first radio peak leads the 7-ray peak by 
phase 0.095 ± 0.002 in phase. Observations of the evolution 
of the peak positions with the energy allow detailed studies of 
the emission regions in the magnetosphere. The delay of the 
radio peaks compared to other wavelengths (optical, X- and 7- 
rays) gives another constraint in the modeling of the emission 
pro cesses taking place in pulsar magnetospheres. In particu- 
lar, [^iperetld] <12QQ3|) reported that, in the fram ework of the 
three-d imensional outer gap model developed bv lCheng et al] 
(2000), one can reproduce the delay of the radio main pulse 
by shifting the production site of the radio emission inwards 
toward the neutron star relative to that of high energy photons. 

In the PC models, 7-rays created near the neutron star 
surface interact with the intense magnetic fields resulting 
in a sharp turnover in the few to 10 GeV energy range, 
while OG and SG models predict a simple exponential cut- 
off. Furthermore, the maximum energy of observed pulsed 
photons must lie below any 7-B pair production turnover 
threshold, providing a lower bound to the altitude of emis- 
sion. We can use the observed phase-averaged cutoff en- 
ergy 6 GeV) to estimate a minimum emission height as 
r > {e max B 12 /l.7SGeV) 2 / 7 P- 1 ' 1 R*, where e max is the 
unabsorbed photon energ y, P is t he spi n period and the sur- 
face field is 10 12 Bi 2 G (Baring 2004). Using the param- 
eters of the Crab Pulsar (P = 33 ms, B 12 = 3.78), 
one obtains r > 3.4i?* which precludes emission near the 
stellar surface. Since we see pulsed photons up to almost 
tmax ~ 20 GeV in the two main pulse peaks, emission at 
these phases must arise at r > 4.8i?*, with a strict lower 
bound of r > 3.7 R* applying to the choice of 8 GeV, the 
lower energy in the highest data point window for the phase- 
resolved spectra in Figure [7] A similar lower bound to the 
emission altitude was recently reported by the MAGIC collab- 
oration using the hyper-ex ponential cutoff e nergy observed on 
the Crab Pulsar spectrum dAliu et al.ll2008l) . We should note 
here that the cut-off energy derived by the MAGIC collabora- 
tion for a simple exponential cut-off (17.7 ± 2.8 ± 5.0) GeV 
is higher than the one obtained with the ferm/-LAT data, 
E c = (5.8 ± 0.5 ± 1.5) GeV. However, the cut-off en- 
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ergy obtained with the LAT using the softer EGRET spectrum 
(7 = 2.022) as done by MAGIC is within the uncertainties of 
the MAGIC value. 

To estimate the pulsed high energy 7-ray efficiency rj 
of a pulsar, one needs to know the total luminosity radi- 
ated L 1 . It can be estimated using L 1 = 4TrffiF }, s D 2 , 
where F obs = (1.31 ± 0.02 ± 0.02) x 10~ 9 erg cm" 2 s" 1 
is the observed phase-averaged energy flux over 100 MeV, 
D = (2.0 ± 0.2) kpc is the distance to the pulsar and /n is a 
correction factor that takes into account the beaming geome- 
try, depending upon the magnetic inclination angle a and the 
Earth viewing angle £ from the rotation axis. For the Crab 
Pulsar, ( is estimated to be (63 ± 2)° from X-ray ob serva- 
tions of the Crab Nebula torus dNg and Romanil 120081) . The 
estimated value of a depends on the e mission model used to 
interpret the data (IWatters et al.1 120091) . The slot gap or two 
pole caustic model best reproduces the observed pulse profiles 
for a ~ 55 - 60°, whereas for the outer gap model a ~ 70° 
gives the best result. Optical polari zation measurements yiel d 
a estimates consistent with these (Slowiko wska et al]|2009l) . 
The corresponding correction factors /n are then equal to 1.1 
for the TPC and 1 .0 for the OG models. This yields a luminos- 
ity of (6.25 ± 0.15 ± 0.15) x 10 35 erg s" 1 above 100 MeV. 
This val ue is consisten t wit h the heuristic lumi nosity law men- 
tioned in lAronsi ( 1 19961) and lWatters et al.l d2009l) . according to 
which r\ oc E~ x / 2 and verified by several 7-ray pulsars such 
as Vela, PSR J202 1 +365 1 (assuming a distance of the order of 
2-4 kpc), Geminga, CTA1, etc. For a neutron star moment of 
inertia of 10 45 g cm 2 , the pulsed high energy 7-ray efficiency 
r) can be derived from the luminosity and the spin down power 
E: j] = L 7 / E = (1.36 ± 0.03 ± 0.03) x 10" 3 above 
100 MeV. 

Knowing the value of the Earth viewing angle £ of (63 ± 
2)°, the main peak separation, which is of the order of 40% 
of the phas e, would be expect ed to be smaller in the Polar 
Cap model (Wa tters et al.ll2.009h . This mismatch persists even 
if moderate altitudes are considered: the open field line cone 
opening angle enlargens to around 12.3° for the Crab at al- 
titudes of around 6 stellar radii, the minimum bound inferred 
above from the observed absence of magnetic pair attenuation 
below 20 GeV This opening angle i s still somewhat too small 
according to the Wat ters et al.1 (120091) analysis to generate the 
observed main peak phase separation. 

6. SUMMARY 

Using 8 months of survey data with the Fermi Large Area 
Telescope and the very precise timing solution provided by 
the Nangay and Jodrell Bank radio telescopes, we have exam- 
ined the high energy behaviour of the Crab Pulsar and Nebula. 
The unprecedented statistics of the data improve our knowl- 
edge of these two sources and place new constraints on theo- 
retical models: 

1 . The 7-ray profile of the pulsar consists of two main peaks 
(PI and P2), very stable in position across the 7-ray energy 
band and separated by S(f) = 0.398 ± 0.003 in phase. The ratio 
P1/P2 decreases with energy, as seen for several other pulsars 
and especially for Vela. The first 7-ray pulse leads the radio 
main pulse by (281 ± 12 ± 21) /zs. 

2. The spectrum of the nebula in the energy range 100 
MeV - 300 GeV is well described by the sum of two power- 



laws of spectral indices T sync = (3.99 ± 0.12 ± 0.08) and 
T IC = (1.64 ± 0.05 ± 0.07) describing the synchrotron 
and inverse Compton components of the Crab Nebula spec- 
trum respectively. No cut-off energy can be estimated for the 
synchrotron component using the LAT data only. The IC ris- 
ing edge studied in the LAT energy range extends nicely up to 
the energy domain covered by Cherenkov experiments. No 
significant cut-off at high energy is observed with the cur- 
rent statistics in the LAT energy range. No significant varia- 
tion in either the synchrotron or Compton components is seen 
with the current statistics on time scales of one, two, or four 
months. 

3. The phase-averaged 7-ray spectrum of the Crab Pulsar 
can be represented by a power-law with an exponential cut-off 
at E c = (5.8 ± 0.5 ± 1.2) GeV. The hyper-exponential cut-off 
index b = (0.89 ± 0.12 ± 0.28) is not significantly favoured 
with respect to the simple exponential b = 1. If only statistical 
errors are included, b = 2 is rejected at 4.9 a level. Using the 
observed cut-off energy to estimate the minimum emission 
height r of the emission region, one obtains r > 3.4i?» which 
precludes emission near the stellar surface. 

4. The pulsar emission is hardest in the phase interval be- 
tween the peaks, usually called the "bridge", while the soft- 
est components is the falling edge of the second peak. Both 
peaks present similar spectral indices, while the cut-off of PI 
is lower than P2, consistent with the energy-dependence of 
the pulse profiles and of the ratio P1/P2. 

5. Knowing the Earth viewing angle £ ~ 63° and the 
value of the inclination angle a comprised between 55 - 60° 
for slot gap models and ~ 70° for the outer gap, one can 
estimate a pulsed high energy 7-ray efficiency of ~ 0.1% for 
the conversion of the spin-down energy to 7-ray emission. 
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FIG. 5. — Spectral energy distribution of the Crab Nebula renormalized to the total phase interval. The tit of the synchrotron (purple dashed line) and IC 
(blue dashed-dotted line) are represented separately with two power-laws. The black curve is the best fit obtained wi th th e sum of these two power-laws. 
The LAT spectral points are obtained using the model-independent maximum likelihood method described in Section l4!2l The statistical errors are shown 
in black, while the red lines take into account both the statistical and systematic errors. Horizontal bars delimit the energy intervals. EGRET data points 
iKuiper et al. 2001) are shown for comparison (green stars). 
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FIG. 6. — Spectral energy distribution of the Crab Pulsar averaged over the whole pulse period. The black curve represents the best fit model, obtained 
with a power-law with an exponential cut-off. The LAT spectral poi nts ( cf. Figure \5\ for the d escription of the conventions) are obtained using the 
model-independent maximum likelihood method described in Section [4~2l EGRET data points (Kuiper et al. 2001) are shown for comparison (green 
stars). 
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FIG. 7. — Phase-resolved spectral energy distributions of the Crab Pulsar. The labels indicate the phase intervals. Spectral results are presented in Tablef2] The 
black dotted curve is the best fit power-law with an exponential cut-off. The LAT spec tral points (cf. Figure[5]for the description of the conventions) are obtained 
using the model-independent maximum likelihood method described in Section B31 90 % C.L. upper limits are computed when the statistical significance is 
lower than 3 <r. 
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FIG. 8. — Variation of the spectral indices, cut-off energies and photon flux above 100 MeV (divided by the phase interval width) as the function of the pulse 
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FIG. 9. — The spectral energy distribution of the Crab Nebula from soft to very high e nergy 7-rays. The fit of the synchrotron component, using COMPTEL 
and LAT data (blue dashed line), is overlaid. The predicted inverse Compton spectra from Atovan and Aharonian 1 1996) are overlaid for three different values of 
the mean magnetic field: 100 fj,G (solid red li ne), 200 fiG ( dashed green line ) and the canonical equipartition field of the Crab N ebula 300 fj,G (dot ted blue line) 
Refere nces: CGR O COMPTEL and EGRET: Kuiper et al. ( 200 1 ^MAGIC: lAlbert et al. (2008); HESS: Aharonian et al. (2006); CANGAROO: Tanimori et al. 
1 1997); VERITAS: Celik (2007); HEGRA: Aharonian et al. (2004); CELESTE: ISmith et all BOO6) 



